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We study the thermal Casimir effect between two thick slabs composed of plane-parallel layers 
of random dielectric materials interacting across an intervening homogeneous dielectric. It is found 
' that the effective interaction at long distances is self averaging and its value is given by a that 

between non-random media with the effective dielectric tensor of the corresponding random media. 
The behavior at short distances becomes random (sample dependent) and is dominated by the local 
values of the dielectric constants proximal to each other across the homogeneous slab. These results 
are extended to the regime of intermediate slab separations by using perturbation theory for weak 
disorder and also by extensive numerical simulations for a number of systems where the dielectric 
function has a log-normal distribution. 



PACS numbers: 03.70.+k, 03.50.De, 05.40.-a, 77.22.-d 



INTRODUCTION 



Systems with spatially inhomogeneous dielectric constants exhibit effective van der Waals interactions arising from 
the interaction between fluctuating dipoles in the system [l], These fluctuation interactions have two distinct 
components: (i) a classical or thermal component due to the zero frequency response of the dipoles and (ii) a 
quantum component due to the non-zero frequency/quantum response of the dipoles. Despite the clear physical 
differences in these contributions, the mathematical computation of the corresponding interaction is almost identical 
and boils down to the computation of an appropriate functional determinant. The full theory taking into account 
both of these component interactions is the celebrated Lifshitz theory of van der Waals - dispersion interactions [H , 
based on boundary conditions imposed on the electromagnetic field at the bounding surfaces and the fluctuation- 
dissipation theorem for the electromagnetic potential operators. From this formulation one can derive the original 
Casimir interaction [4J by taking the limit of zero temperature and ideally polarizable bounding surfaces. In this 
respect the Lifshitz theory is nothing but a proper finite temperature and realistic boundary conditions generalization 
of the Casimir interaction. The Lifshitz - van der Waals interactions is thus indeed nothing else but the thermal 
' Casimir effect. 

The major mathematical problems in the computation of Casimir type interactions (setting aside the experimental 
and theoretical challenges to determine the correct dielectric behavior) are (i) the application of the Lifshitz approach 
to non-trivial geometries (i.e. beyond the cases of planar, spherical and cylindrical geometries) and (ii) taking into 
account local inhomogeneities in the dielectric properties of the media, always present in realistic systems and thus 
relevant for the comparison of theory with experiment. In this paper we will address the second of these points in 
detail, and to our knowledge we present the first analysis of the effect of dielectric disorder on Lifshitz - van de Waals 
interactions (apart from a recent letter [j| by the authors on this subject). Specifically, we will consider the thermal 
(zero-frequency) Casimir interaction for the case where the local dielectric constant is a random variable. Specifically 
we will consider the interaction between two thick parallel dielectric slabs, separated by a homogenous dielectric 
medium, see Fig. (JTJ. The dielectric response within the two slabs is constant in the planes perpendicular to the slab 
normal, but varies in the direction of the surface normal. It is well known that this problem can be solved in the 
case where the dielectric constants of the slabs do not vary [2| and the result can be tentatively applied to the case of 
fluctuating dielectrics constants via an effective medium theory which consists of replacing the fluctuating dielectric 
constant by an effective (spatially homogeneous within each of the slabs) dielectric tensor. Naively one might try the 
approximation 

e(x)-<e>, (1) 



FIG. 1: (Color online) A schematic presentation of the model. Two finite slabs with disordered plane-parallel dielectric layers 
interacting across a dielectrically homogeneous slab of thickness £. z axis is perpendicular to the plane of the slabs. 



where the angle bracket denotes the spatial or ensemble averaged dielectric constant within the slab in question. 
However the most commonly used approximation is that, where the local dielectric tensor is replaced by the effective 
dielectric tensor P, 0], i.e. 



where the bulk dielectric tensor is defined via 



e«(x)-4 e) > ( 2 ) 



e^{E j )={e ij E j ). (3) 



The use of the effective dielectric constant is not easily justifiable mathematically as an approximation, although 
physically the effective dielectric constant clearly does capture the bulk response to constant electric fields. We 
shall see in this paper, for the random layered dielectric model studied here, that the effective dielectric constant 
approximation of Eq. ^ does in fact give the correct value of the thermal Casimir interaction when the two slabs 
are widely separated. On can argue that this is to be expected on physical grounds as the fluctuating electromagnetic 
field modes with small wave-vector (corresponding to variations on large scales) dominate the Casimir interaction 
for large inter-slab separation. The dielectric response of the material to a constant electric field is given by the 
effective dielectric constant and if the wave-vector dependent response is suitably analytic near k = we expect that 
eg>(k)~ e W(0) = e«for|k|«l. 

II. THE MODEL AND GENERAL ANALYSIS 
A. Formulation 

The Hamiltonian associated with the thermal fluctuations of the electrostatic field in a dielectric medium is given 
by the classical electromagnetic field energy 

H\d>] = \ J dx e(x) (V^x)) 2 (4) 
and the corresponding partition function is given by the functional integral 

Z = J <#]exp(-/3#[$). (5) 
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Differences in dielectric constants lead to the thermal Casimir effect which arises from the full treatment of the 
thermal (zero frequency) van de Waals forces in the system. Here we will consider layered systems where the dielectric 
constant e only depends on the z direction e(x) = e(z). If we express the field <j) in terms of its Fourier modes in 
the plane perpendicular to z, coordinates denoted by xi and which we will take to be of area A, with wave- vector 
k = (k x , ky) then the Hamiltonian can be written as 



(6) 



with 



#k = - I dz e(z) 



^(z,k)^(z,-k) +k2 ^ kfc _ k) 



dz 



dz 



A direct consequence of this decomposition of the Hamiltonian is that the partition function can be expressed as a 
sum over the partition function of the individual modes as 



ln(Z)=]Tln(Z k ) 



(8) 



where 



Z k = J d[X] exp ^-i J dz e(z) 



^] + k 2 X 2 
dz 



(9) 



Here k = Ikl and we have taken into account that the field 6 is real. 



B. Evaluation of the functional integral 



The problem of computing the interaction between slabs composed of layers of finite thickness can be studied using 
a transfer matrix like method Q. However we will use a method based on the Feynman path integral instead, which 
is particularly well suited to the study of systems where the dielectric constant can vary continuously in only one 
direction [jj • If we specify the starting and finishing points of the above path integral to be x and y respectively at 
times z' and z we see that it has to be of harmonic oscillator form defined by 



K(x,y;z',z) 



X(z)=y 



d[X] exp 



X(z')=x 



dz M (z) 



dX\- 
dz J 



L0 2 X 2 



(10) 



where the mass, which is z dependent, is given by M(z) — e(z) and the frequency lo is given by u> = k. In the case 
where M and uj are constant, the propagator K is given by the well known formula 



K(x,y;z',z,M,uj) = 



Muj 



2-7T sinh(w(z — z')) 



exp 



- ((x 2 + y 2 )Mujcoth(u(z - z'))2xyMujcosech(uj(z - z'))) 



In the case where M (or indeed u>) vary with z we can still formally compute the path integral via the generalized 
Pauli - van Vleck formula which tells us that K must have the general form 



K(x,y;z',z) = [ — 



exp 



~a,i(z', z)x 2 - ^a f (z', z)y 2 + b(z' , z)xy 



(12) 



We may now write down an evolution equation for the coefficients Oj, aj and b using the Markov property of the path 
integral (in fact this is how one can prove the generalized Pauli - van Vleck formula (7[) 



K(x, y; z' , z + () — / dw K[x, w; z' , z)K(w, y; z, z + Q. 



(13) 



Now if we take £ = Az infinitesimal and assume that M and u> are constant over (z, z + Az) (but could have jumped 
at z), then by looking at the coefficients of x 2 , y 2 and xy in the so computed K(x, y; z + Az) we find the following 
evolution equations for a^, a/ and b : 
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da.i(z',z) b 2 (z',z) 

Oz M ' K ' 

db(z',z) b(z',z)a f (z',z) 

~8z~ ~ M ' (15j 

Note that the evolution equation for b can also be obtained by examining the change in the pre-factor term of the 
propagator K. An equivalent and related way of deriving these equations is to note that, from the Feynman - Kac 
formula, K(x, y; z, z') satisfies the Euclidean version of the Schrodinger equation: 

^^'^K^^-^H^^' 4 (17) 

Substitution for K(x, y; z' , z) from Eq. (fP2|) then yields the evolution equations ((T6|) . As the action is positive definite 
we expect that both m and a/ are positive (by considering paths staring or ending at 0). Also for \z — z'\ large, 
with respect to the correlation length of the disorder, we expect K to factorize in its x and y dependence and thus 
the coefficient b should decay to zero for sufficiently thick slabs. Note that if e is a stationary process then the path 
integral should be the same run backwards in time, z, as when it is run forwards. This means that ai and aj should 
have the same steady state distribution. 

Consider a system of two thick slabs of respective thicknesses L\ and Li separated by a distance I and the region 
between them occupied by a dielectric medium of dielectric constant eo - a vacuum or air, for example. From our 
discussion above, for large L\ and Li the partition function for the mode k is thus proportional to 



— — ((.t 2 + y 2 )eok coth(fcz) — 2x?/eofccosech(fcz)) 



x exp(-i4 2) (fc)y 2 ), (18) 

where (k) is the solution of Eq. (116|) with some initial conditions (which for an infinite slab do not have to be 

(2) 

specified) given at z = z = —L\ to zero with to — k and M{z) = e(z) and aj (k) is the corresponding quantity for 
the slab 2 (with initial conditions given at z = z' = Li + I and evaluated at z = I). 

We thus find that the ^-dependent part of the free energy of the mode k (up to a bulk term which can be subtracted 
off to get the effective interaction.) is given by 

„ k BT] ( {a { P(k)-e k){af{k)-e Q k) \ 

F k = — - In 1 J — J m exp -2fcZ , (19) 

2 V {af{k) + e Q k){af (k) + e k) J 

with the total I dependent free energy given by 

F = VF k . (20) 



In order to evaluate the integrals of a} (k), one first has to solve equations of motion Eqs. (TT6|) to get the z 
dependence of af(k,z) and then proceed to the integrals that enter Eq. (fl9|) . The evolution equation for a/(fc) for 
either slab can be read of from Eq. (fT!))) and is given by 

dz e(z) 

where we have dropped the explicit dependence on the initial point z' — —L%. Though it does not look like it at the 
first sight, this equation is simply a rewriting of the underlying Poisson equation for the original, charge free, dielectric 
system . This can be seen as follows. Assume first of all that a/(/c, z) can be parameterized with a function ip(k, z) as 

a f (k,z) = e(z)^lnip(k,z). (22) 
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In quantum mechanics problems with disorder the above change of variables is often used since in the presence of 
disorder nonlinear first order equation is easier to analyze then the second order linear equation [8[. Inserting now 
this ansatz back into Eq. ([2Tj) we find that it implies 

5 0£) < 23 > 

for ip(k,z), which, is nothing but the Poisson equation for an inhomogeneous dielectric, where the inhomogeneity is 
only in the z direction, which has been Fourier transformed in the directions perpendicular to z. The a/(fc,z) in 
Eq. (|22|) is thus given by the solution of the Poisson equation in the specified planar geometry. This is of course no 
surprise since we are indeed dealing with an inhomogeneous electrostatic problem. On the other hand the derivation 
presented above is completely equivalent to the transfer matrix method [6j or to the density functional method [9j for 
evaluating the van der Waals forces. One of the clear strengths of this method is that it allows the of computation 
the van der Waals interaction to be carried out using a local method where the coefficients cif, cif and b for any of the 
media involved can be computed and then the interactions between any combinations of media can be worked out in 
terms of these coefficients. 

If we now write ap(k, z) — ka^(k, z) and if the distributions of the a^(k, z) — y are given by Pi(k, y) then we 
find that, in three dimensions the average I dependent free energy is given by 

(F) = JdkkJd Vl J dy 2Pl (k,y 1 )p 2 (k,y 2 )ln (l - & ~ ^ ~ ^ exp(-2fc/)) , (24) 

where the angled bracket on the l.h.s. indicates the disorder average over the dielectric constant within the slabs and 
we have assumed that the realizations of the disorder in the two slabs are independent. This is why the joint disorder 
probability distribution is multiplicative for the two layers. 

In the case where the slab thickness is large, L — ► oo, it is simple to derive a scaling formula from Eq. (j2Tj) for the 
probability distribution. Let e(z) = e g f(z) where f(z) is an instance drawn from an ensemble distributed according 
to a given distribution. Then p{k 1 y\ e g ) is the resulting distribution for y given the mode wave-vector k. We then find 
that 

p(k, y; Xe g ) = jP( k , v/\ e 9 ) • (25) 

Thus, we need to compute the probability distribution once only, say for e g — 1, and obtain those for other values of 
e g using this scaling result. 



C. Large I limit 

Let us first investigate the form of van der Waals interaction free energy in the limit of large separations between 
the two slabs. The equation obeyed by a is 

da(k,z) ka 2 

— j = e« , (26) 

dz e 

which can be written as 

d u( k , ,,,,, Q 2 , 07 ^ 

= e{c/k) ~ WJkY (27) 

with £ = zk. When k is small e(£/fc) varies very rapidly and thus becomes de-correlated from the value of a. The 
Laplace transform for the probability density function of a is defined by 



p(k,s,()= dyexp(-sy)p(k,y,C) = (exp(-sa(fc, £))) , (28) 
Jo 

and, from the equation of motion Eq. (I26|) . obeys 

J ' /i " / '" S ° ' d :(exp(- S «(fc,C))> =/e(CA)exp(- S a(fc,C))--^ 7T exp(- S a(fc,C))V (29) 
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Assuming that k is small and thus a(k, £) and e(^/k) are de-correlated we can write 

- = - (Ve) ^P(k,s,0 • (30) 

This equation has the large £ equilibrium solution (justified as we are taking the limit L\, L2 — > 00) 

lim p(k, s, C) = exp(-e*s) , (31) 

c — >0 ° 

with 



Inverting the Laplace transform then gives the equilibrium distribution 

p(y,k)=S(y-e*) (33) 

at small k. When I is large the integral in Eq. (|24|) is dominated by the small k behavior and we may use the analysis 
presented above, especially Eq. (|33p in Eq. (|24p. to give the following asymptotic form for the interaction free energy 



F !/ _ X } ~ rS^ l udu ln(1 " A I A 2 ex P(-")) ~ "^T > l84 > 



with 



A* = ^— ^ - (35) 

and where e* are defined via Eq. (|3"2|) . The subscript i on the angled brackets signifies that we are averaging the 
dielectric constant in the slab i. The term H* defines an effective disorder-dependent Hamakcr coefficient. 

This result can be obtained via another more physical argument. The random layered materials can be replaced 
by an effective anisotropic medium where the dielectric tensor is non isotropic and has the form 



eg = e,| , (36) 

e xx = £ yy — e i i (37) 

all other terms being zero by symmetry. The term e| is the dielectric constant in the z direction is known exactly 
and is given by 

,(<0 _ 1 



C = 7T7T > (38) 



(Ve) ' 

and the perpendicular component is simply given by 

ef = ( £ ) . (39) 

The forms of el and ej 6 follow simply from the fact that in the perpendicular direction the dielectric constant is 
obtained by analogy to capacitors in series and in the parallel direction by analogy to capacitors in parallel arrangement 

It is a straightforward exercise to see that the effective value of e* for this system coincides with that of Eq. (|3^|) 
above. This correspondence makes physical sense since the high wave length or small k fluctuations of the electric 
field are responsible for the behavior of the Casimir interaction at large distances and the effective dielectric response 
at low (but non-zero k) must be close to that of the response to a constant field i.e. described in terms of a dielectric 
constant. An interesting consequence of this result is that for large separations (where I is much greater than the 
correlation length of the dielectric disorder) the thermal Casimir interaction free energy is self averaging. 

It is instructive to compare the result for the Casimir interaction at large separations in the case of a fluctuating 
dielectric constant and in the case of a homogeneous medium, whose dielectric constant is the average of that in 
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the fluctuating medium: = (e) (the subscript h signifying that the medium is homogeneous). The homogeneous 
medium has a Hamakcr constant 



16 

with 



— J u du In (l - Al cxp(-u)) , (40) 



Ah W^~o> (41) 
when the medium (1) has the same composition as the medium (2). Jensen's inequality tells us that 



> A (42) 



since the function f(x) = 1/x is convex. Thus 



e * = V<w- (e) (43) 

Clearly the effective Hamaker constant is a monotonic function of A 2 and the interaction is always attractive. The 
difference in A 2 for the two systems is 

A *2 A 2_ 1 eo( e *-(e))(e*(e)-e 2 ) 

h ( ( £ ) +eo) 2 (e * +eo) 2 ■ ^ > 

Therefore, using Eq. (|43|) . we find that, if e*(e) > eg then > H* , and the interaction between the two 
homogeneous media is stronger than that between the two fluctuating media. This condition can be written as 

' ' >4 (45) 



and is always satisfied if (l/e) -1 > eo. However it is always violated (the interaction between the homogeneous 
media is weaker than that between the random media) if (e) < eo- We thus see that, depending on the details of 
the distribution of the fluctuating dielectric response in the two slabs and the dielectric response of the medium 
in-between, the effective interaction at large inter-slab separations can be stronger or weaker than that for a uniform 
medium with a dielectric constant equal to the average value of that of the media in which it fluctuates. 

D. Small I limit 

One would imagine that as the distance between the slabs is reduced the result will be increasingly dominated by 
the slab composition at the two opposite faces 0]. Indeed in the small I limit Eq. ([24]) is dominated by the large k 
behavior. The asymptotic behavior can be extracted if one assumes the ansatz 

«(M) = E^- (46) 

n=0 

Substituting this into Eq. (f2"6")) gives the following chain of equations for a n (z) 



1 1 da n (z) _ 1_ y, a n (z)a m (z) 

k^ k n dz ~ £[Z) e(z) ^ k m + n ' 1 ' 

n— n,m— 

^From here it is easy to see that to order O(l) the leading asymptotic result of Eq. ([52]) is given by 



a (z) = e(z) 



(48) 
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The equation for the corrections (n > 1) to this asymptotic limit is 

d(Xn-l(z) _ 1 



dz e(z) 

v ; m=0 



ot m {z)a n - m {z), (49) 



and the next two terms in the asymptotic expansion are given by 



4 dz 2 8e(z) \ dz 

It is straightforward to realize that these terms generate 0(1/7) corrections to the asymptotic result which are 
subdominant when / is large. Thus to the leading order 

a(z, k) « oto(z) = e(z) (52) 

and from here it follows straightforwardly that 

lim Pi(y,k) = p % (e) (53) 

k — >oo 

where pi is the probability density function of e(z) in medium i. This result is easily understood from the physical 
discussion above. The average of the thermal Casimir interaction free energy Eq. (|24[) in the small separation limit 
is thus given by 

(F)(Z— +0)~ J udu J de 1 de 2Pl (e 1 ) P2 (e 2 ) ln(l- A X A 2 exp(-u)) , (54) 

with 

A s: = ^^. (55) 
e* - eo 

The forms of the thermal Casimir interaction free energy are thus given by Eqs. (|54[) and (|34|) in the small and 
large inter-slab separation limits respectively. We have thus obtained the limiting behavior of the thermal Casimir 
effect in the limit of large separation between the slabs, where the free energy is given by self-averaging and thus the 
distributions of a(k,z) are strongly peaked, and in the limit of small separation, where the free energy is a random 
variable that has to be averaged over the probability density function of the dielectric constant in the media composing 
the two interacting slabs. 



III. SMALL DISORDER LIMIT - PERTURBATION THEORY 



The analysis presented above is valid for any type of disorder, irrespective of its properties. Here we also investigate 
a different approach that takes into account the disorder effects on a perturbative level i.e. another way of approaching 
this problem analytically is to assume that the disorder is small. We assume in this case that the dielectric response 
can be written with an ansatz of the form 



e(z) = e g exp(X(j>(z)) , (56) 

where A is a scalar parameterizing the strength of disorder which will be used as the expansion parameter for per- 
turbation theory. When A = we have the disorder free homogeneous system which is the starting point for the 
perturbation expansion (zeroth order). If the mean of the field of <f> is zero then e g is the geometric mean of the 
dielectric constant. We now assume that A is small and write 



oc 

a(z, k) = X n a n (z, k). 

?i=0 



(57) 
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Substituting this into Eq. ((26)) and matching the powers of A, we obtain the following equations for the first three 
terms in the perturbation expansion when z is large: 



fc(-2ai + 2e g cj>) , 



a = €g , 

da i 

dz 

da2 , , „ a? „ ,s 

— - = fc(-2ar 2 - — + 2ai$. 

dZ €g 

At large z we can solve the last two equations to obtain 

a\ = 2kt g I dz' exp(— 2k(z — z'))4>(z') , 
Jo 

a 2 = 4k 2 e g [ dz' exp(-2k(z - z'))<l>(z') [ dz" exp(-2k(z' - z"))<j>{z") 
Jo Jo 

z r ' i 2 

- 4k 3 e g [ dzexp(-2k(z- z')) [ dz" exp(-2k(z' - z"))<f){z") 
Jo Jo 



(58) 
(59) 

(60) 



(61) 



(62) 



We can now verify some of our previous results at the level of perturbation theory. For k — > oo we have that 
2kexp(— 2k(z — z')) — » 8{z — z'), and this gives 



This gives 



ai -» e g cj)(z) , 
aa - 1 2 W. 



a e s (l + A0 + y</« 2 ) w e g exp(A0) = e(z) 



in agreement with the large k result stated previously. Also in the limit of small k we will have that 



lim 2k / dz' exp(-2fc(z - z'))4>{z') -> («/>) , 



for large z and so we have 



Since (<f) =0, we find 



ai -> e s (0) , 

«2 - |(0> 2 - 



lim a(z, k) = e g , 

k^O 



which is in agreement (to 0(A 2 )) with the general result Eq. 

Now let us consider the case where the field <fr is Gaussian of zero mean with correlation function 

{4>{z)4>V)) =exp(-cj\z-z'\), 



(63) 
(64) 

(65) 

(66) 

(67) 
(68) 

(69) 
(70) 



where 1/oj is the correlation length of the field. This Gaussian field has the correlation function of an Orstein- 
Uhlenbeck (OU) process and is Markovian. The moments relevant to the degree in perturbation theory we are working 
to (0(A 2 )) are 



(ai) = 0, 



«i = = 



(a 2 ) 



2k + uj 
2k + u> ' 



(71) 
(72) 

(73) 
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FIG. 2: Probability density for a for k = 1, 0.1, 0.01 for e(z) = exp((f>(z)) where <j> is the OU Gaussian process with correlation 
function given by Eq. (|70[) . The distribution p(y,k) becomes sharply peaked around y = 1 which is as predicted for small k, 
e* given by Eq. J32} 



Using this we may write the corresponding free energy as a random variable 

k B T 
167rZ 2 



F = TF~?2 I udu]n(l- A[A' 2 exp(-u)) , ( 74 1 



where 



with 



4 - Jr-T 1 , (75) 

£j + £0 



~9,» 



i + V— JL -7^ + ^~ ^ ) ■ ( 76 ) 
w u + ujil 2 u + 0Jil 1 



The subscript i refers to the medium (1, or 2) and the <7j are independent Gaussian random variables of zero mean 
and unit variance. Note that to order A 2 we can replace the term «2 by its mean and have done so in the above. 
Also we should expand the corresponding expression for the free energy to order A 2 , this will ensure that the resulting 
average and moments will always be finite. 



IV. NUMERICAL SIMULATIONS 



In this section we perform some numerical simulations to verify the general asymptotic behavior of the Casimir free 
energy that we have analyzed in earlier sections. For large slab separation (large X) the results are summarized by 
Eqs. ([55)1 and (f52")) . and in the small separation limit (small I) we expect Eq. ([55)1 to hold. We present simulations to 
show that these equations are indeed correct in the relevant limit. 

We recall that for small k the values of evi and «2 are predicted to be self averaging and thus their distributions 
should be strongly peaked as k — > 0. 

In Fig ^ we have simulated the system when e is given by Eq. (|5"6"|) and 4> is an OU Gaussian process (with A = 1 
and e g = 1 in Eq. (|56|) and id = 1 in Eq. ([70)) ) for values of k = 1 to 0.01. We see that as predicted the distribution 
of a becomes increasingly peaked about e* = 1 (as given by Eq. ([32]) ) as the value of k is decreased. For large k the 
prediction Eq. (|53|) can be verified. Shown in Fig. (j4]) is the distribution for the same distribution of e(z) but for 
k = 10, we see that, as predicted, the probability density function for a is already very close to that of e which is 
given by 
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2.5 




FIG. 3: Probability density for a for k = 1, 0.1, 0.01 for e(z) — exp(<f)(z)) where (f> is the Gaussian process with correlation 
function given by Eq. (|78[) . As for the OU process the distribution p(y,k) becomes sharply peaked around y = 1 which is as 
predicted for small k, and e* is given by Eq. (1321) 



Note that q{y) is independent of the correlation function of the Gaussian field in the log-normal distribution. 

One can also consider the case of non-Markovian log-normal dielectric constants. For instance one can take the 
field </> in Eq. (|56[) to have correlation function 

(0(z)0(z'))=exp(^ 2 (z-z') 2 /2). (78) 

In what follows we fix A = 1 and to = 1. For the small values of k shown in Fig. ([3]) we see again that as k — > 
that the distribution becomes peaked about e* = 1 as predicted. In Fig ([5]) we show the comparison between the 
distribution of a and e, again at this value of k the agreement between the two distributions is already excellent. 

We can also compute the average effective Hamaker coefficient as a function of separation I using Eq. (fM| to define 

(H(l)) ^^-JuduJ dy 1 dy 2Pl (u/2l,y 1 )p 2 (u/2l,y 2 )]M(l-A(y 1 )A(y 2 )exp(-u)) , (79) 

with A(y) = (y — £o)/(y + £o)- For e(z) distributed according to the log-normal distribution in Eq. 156|) and using 
the scaling result in Eq. (|25[) we can write 



(H(l;e g )) = J udu J dy x dy 2 pi{u/2l, y{)p 2 (u/2l, y 2 ) m(l - A(e g yi) A(e g y 2 ) exp(-u)) 



(80) 



where the Pi(k, y) are computed for e g — 1. We show H(l; e g ) versus I in Figs. ((6]) and (|7|) for various values of e g . 

As can be seen H(l) is asymptotic to the Hamaker constant for a homogeneous system with dielectric constant 
e* = yj (e)/(l/e), and that for e g = 1.0 that H(l) is asymptotic to zero from above showing that the force is attractive 
for all I. This need not have been the case since there will be contributions from configurations were the dielectric 
constants ej., e 2 in the slabs satisfy t\ > 1 > e 2 or vice- versa; such configurations contribute a repulsive contribution to 
the force. As is also observed, the curves are either asymptotic from above or below and are not necessarily monotonic 
but can have an initial decrease before rising to the asymptotic value. These are new and significant features due to 
the random nature of the system. 

In Figs © and © we show a sample of the Hamaker coefficients for particular realizations of e(z) plotted against 
separation I and for e g = 2.0, 10.0, respectively. The averages over the ensembles from which these curves are taken 
are given in Fig. (|7|). It is important to note that the curvatures of the ensemble curves are not of definite sign and 
that for small I there are both curves that decrease and curves that increase with /. These properties are reflected in 
the shape of the curve for the ensemble average. 



V. CONCLUSIONS 



We have obtained the limiting behavior of the thermal Casimir effect in the limit of large separation between the 
slabs, where the free energy is given by self-averaging and thus the distributions of a(fc, z) are strongly peaked. We 
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FIG. 4: The filled circles are the probability density for a for k = 10 for e(z) = exp(^>(z)) here <j> is the OU Gaussian process 
with correlation function given by Eq. (|70|l with w = 1. The solid line is the probability density function q(y) given in Eq. (177(1 
for the same distribution of e. The two distributions are already very close for k = 10 in agreement with the prediction of Eq. 




FIG. 5: The filled circles are the probability density for a for k = 100 for e(z) = exp(0(z)) here cf> is the Gaussian process with 
correlation function given by Eq. (|78p with w = 1. The solid line is the probability density function q(y) given in Eq. (|77p 
for the same distribution of e. The two distributions are already very close for k = 10 in agreement with the prediction of Eq. 
Q53). 

have shown that the interaction between two homogeneous media is stronger than that between the two fluctuating 
media if (l/e) -1 > e . However it is always weaker if (e) < e . We thus see that, depending on the details of 
the distribution of the fluctuating dielectric response in the two slabs and the dielectric response of the medium 
in-between, the effective interaction at large inter-slab separations can be stronger or weaker than that for a uniform 
medium with a dielectric function equal to the average value of that of the media in which it fluctuates. In the limit of 
small separation, where the interaction free energy is a random variable that has to be averaged over the probability 
density function of the dielectric functions in the media composing the two interacting slabs. 

The intermediate length scales were analyzed via perturbation theory and models of disorder that can be treated 
numerically. All numerical simulations completely corroborate the analytical results for self-averaging at large sepa- 
rations 

The non-linear formulation of the problem presented here should be equally useful to treat the case of deterministi- 
cally varying dielectric functions and could open up a useful computational framework for designing materials where 
the effective interaction can be tuned, to induce attractive or repulsive forces depending on the separation [111 ]. The 
formulation also means that if one knows the coefficients aj(fc), a/(fc) and bj(k) for any set of slab media, then one 
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FIG. 6: The effective Hamaker coefficient H(l) defined in Eq. (|80[) for the distribution defined by Eqns. (|56[) and (|78|l for 
e g = 1.0, 1.414. In both cases the curves are asymptotic to the Hamaker constant for a homogeneous system with dielectric 
constant 

e * = V(e)/(l/e). Notably, e* = 1 for e g = 1.0 and we see that H(l) is asymptotic to zero from above showing that 
the force is attractive for all I even in this case. As is also seen, the curves are asymptotic either above or below depending on 
the value of e 9 and is not necessarily monotonic. 
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FIG. 7: The effective Hamaker coefficient (H(l)) defined in Eq. (l80|l for the distribution defined by Eqns. ([56} and (f78|) for 
e fl = 2.0, 10.0. In both cases the curves are asymptotic to the Hamaker constant for a homogeneous system with dielectric 



constant € 



= y/(e)/(l/e). 



can immediately compute the effective interaction between them at any distance. This is a rather surprising result as 
if one wanted to compute the effective interaction between two media (1) and (2) using the pair wise approximation, 
1/r 6 for the van der Waals interaction, it is clear that medium (1) needs to know what is in medium (2) at each point 
in order to compute the force. The decomposition in terms of Fourier modes however means that the interaction 
between the two media is effectively factorized. A number of further applications of our method would be to examine 
the role of disorder for the non-zero frequencies corresponding to quantum fluctuations and also for different geome- 
tries, such as cylindrical and spherical, when the dielectric function varies only in the radial direction. An interesting 
open problem concerns what happens when the dielectric function varies in all directions, it would be interesting to 
know if one can prove in this case whether the long-distance interaction between two slabs is also given by the same 
effective medium expression as derived here. 
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FIG. 8: Examples of the Hamaker coefficient from the ensemble generated by the distribution for e(z) for e g — 2.0 plotted 
versus separation I. The average curve of the ensemble is shown in Fig. ([7]). Note that for small I the curvature of the ensemble 
curves are not of definite sign 
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FIG. 9: Examples of the Hamaker coefficient from the ensemble generated by the distribution for e(z) for e g = 10.0 plotted 
versus separation I. The average curve of the ensemble is shown in Fig. @. Note that the curvatures of the ensemble curves 
are not of definite sign. 
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